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Abstract. The Casimir effect for general Robin conditions on the surface of a cylinder in D- 
spacetime dimensions is studied for massive scalar field with general curvature coupling. The 
energy distribution and vacuum stress are investigated. We separate volumic and superficial 
energy contributions, for both interior and exterior space regions. The possibility that some 
^ ■ special conditions may be energetically singled out is indicated. 

1 Introduction 

t-H ! 

The whole scientific production related to the Casimir effect has already a quite impressive 
volume (for reviews see e.g. refs.Jl]). This effect describes the energy variation undergone by 
a quantum system in which the field modes are constrained by specific conditions on a given 
boundary. The most often studied types of boundary and conditions are those associated to well 
known problems, e.g., plates, spheres, and vanishing conditions, perfectly conducting conditions, 
etc. 

| Our aim in the present paper is to analyze a combination of boundary geometry and condi- 

tions which hasn't attracted the researcher's attention up to now. We are referring to general 
Oh. homogeneous — or 'Robin' — boundary conditions on infinitely long cylindrical surfaces. Such 

conditions are an extension of the ones imposed on perfectly conducting boundaries and may, 
in some geometries, be useful for depicting the presence of interfaces between dielectric media. 
Robin type of boundary conditions also appear in considerations of the vacuum effects for a 
confined charged scalar field in external fields || and in quantum gravity ||]. Unlike 'purely- 
Neumann' boundary conditions, they can be made conformally invariant. The Casimir energy 
for cylindrical geometry can be important to the flux tube model of confinement [||, |6| and for 
determining the structure of the vacuum state in interacting field theories Jj7j . The cylindrical 
problem with perfectly conducting conditions was first considered in ref.|8|| (see also refs.||, 
[|10|]). While the earliest studies have focused on global quantities, such as the total energy 
and stress on a shell, the local characteristics of the corresponding electromagnetic vacuum are 
considered in [11] for the interior and exterior regions of a conducting cylindrical shell, and in 



X 

5-h 



[p.2[| for the region between two coaxial shells (see also []13| ), 

In this paper we will study the vacuum expectation values of the energy-momentum tensor 
for the massive scalar field with general curvature coupling satisfying Robin boundary condition 
on the surface of a cylindrical shell in arbitrary spatial dimensions. As for other geometries, an 
important question (see, for instance, |l4j]) is the relation between the mode sum energy, as a 
renormalised sum of the zero-point energies for each normal mode, and the volume integral of 
the renormalised energy density. We will show that for the geometry under consideration these 
two quantities are different and will interpret this difference as a result of the additional surface 
energy contribution (similar considerations for the plane and spherical geometries were made in 
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The paper is organized as follows. Sec. 2 is devoted to the energy density inside a cylindrical 
shell, from the expectation values of the energy-momentum tensor. Calculations of this nature 
involve infinite sums. The summation method applied in this case is based on a variant of the 
generalized Abel-Plana formula [|l^] together with an adequate cutoff function. In sec. 3, the 
integrated energy inside the cylinder is calculated. One finds that, in general, the total vacuum 
energy differs from the integral of the volumic density. This is explained by the presence of a 
surface contribution located at the boundary. The part corresponding to the modes propagating 
in the exterior is considered in sec. 4. Then, it is possible to sum both pieces and obtain a 
formal expression for the total integrated Casimir energy. In sec. 5 we study the possibilities 
of performing actual evaluations of this magnitude. Zeta function is proposed as a summation 
method, and its working is illustrated in the example of D = 3 + 1 and zero mass. The ensuing 
ending remarks are included into sec. 6 



2 Casimir energy density and stress inside a cylindrical shell 

We will consider a scalar field (p, with curvature coupling £, satisfying a generic Robin boundary 
condition 

(A + BrtV^tpix) = (2.1) 

on the cylindrical shell with radius a. Here n l is the normal to the boundary surface, V« - 
is the covariant derivative operator, A and B are constants. The results in the following will 
depend on the ratio of these coefficients only. However, to keep the transition to the Dirichlet 
and Neumann cases transparent we will use the form (|2.1| ). The corresponding field equation is 
in form 

(V*V< + £R + m 2 ) tp[x) = 0, (2.2) 

where R - is the scalar curvature for the background spacetime, m - is the mass for the field 
quanta. By using this equation, in the case of the flat background the corresponding metric 
energy-momentum tensor (EMT) may be presented in the form 



(£ " ~)9ik v'Vj - £V,;V fc 



^. (2.3) 



The vacuum expectation values (v.e.v.'s) for these quantities can be derived by evaluating the 
mode sum 

(0|r ifc (x)|0) =Y,T ik {<p a {x),<p* a (x)}. (2.4) 

a 

Here {(p a (x)} is a complete orthonormal set of positive frequency solutions to the field equation 
with quantum numbers a, satisfying the boundary condition (2.1). 

In this section we will consider the scalar vacuum inside a cylindrical shell in D-dimensional 
spacetime. In accordance with the problem symmetry we will take cylindrical coordinates 
(r,4>,Zi,..., zn), N = D — 3. The corresponding eigenfunctions to ( |2.2| ) have the form 



<-Pa{x) = Pa Jn {ir) exp (in(j) + zkr || — iiot^j , uj = \J"f 2 + k 2 + m 2 , (2-5) 
a = (n, 7 ,k), n = ..., -2,-1, 0,1, 2, • • • , (2.6) 

where rn = (zi, . . . , zn) an d the coefficients (3 a are determined from the orthonormality condi- 
tion. The eigenvalues for the quantum number 7 are quantized by the boundary condition fl2.1| ) 
on the cylinder surface r = a. From this condition it follows that the possible values of 7 are 
solutions to the equation 

J n ( 7 a) = AJ n {na) + S 7 j;( 7 a) = 0. (2.7) 



2 



Here and below we use the notation 

f(z) = Af(z) + (B/a)zf'(z) 



(2i 



for a given function f(z). It is well known that for real A and B the all zeros for J n {z) are simple 
and real. Let z = \ n j, I = 1, 2, ... be the corresponding positive zeros, arranged in ascending 
order: A n> ; < A ni / + i. Now the possible values for 7 can be expressed as 7 = \ rh i/a. 
From the orthonormality condition 

dVip a (x)tf* a ,(x) = —8 nn >8i V 5{k - k'), (2.9) 
L<jj 

where the integration goes over the region inside a cylinder, one finds for the coefficients (3 a 

,2 _ l T n(ia) 



1% 



ua{2TT) N + 1 ' 



and, following [ 13 1 , we have introduced the notation 

T n {z) 



(z*-n*)J*(z) + z*J<?{zy 



(2.10) 



(2.11) 



Using eigenfunctions (|2.5| ) for the v.e.v. of the two-field product (positive frequency Wightmann 
function) we have 



(0\<p(x)<p(x')\0) = 52<p a (.x)<P*a&) 



1 



E 



(^) N+1 a 2 V ^j\ 2 Ja 2 + k 2 + m 2 
x Jn(K,ir/a)J n (X nj ir'/a) exp \in{(j) - (jJ) + ik(r|| — r'j) — iu(t - t') 



(2.12) 



where 



+00 00 



E = K k E E- 



n=— 00 i=l 



Substituting this expression into (|2.4D for the v.e.v. of the EMT one finds 

(0|2*|0) = diag(e, -p u -p 2 , -p 3 , ■ -PD-i)- 



(2.13) 



(2.14) 



Here e is the vacuum energy density, pi, p 2 , P3 = Pi = ■ ■ ■ = Pd-i are effective pressures in the 
radial, azimuthal and longitudinal directions, respectively (vacuum stresses). These quantities 
are determined by the relations 



l(r) 



1 



E 



fn\ J n{K,ir / a)), q = e,pi 



( 2 ^) JV+1 « 3 a ^j\ 2 n j + k 2 a 2 + m 2 a 2 
where, for a given function f(z), we have introduced the notations 



(2.15) 



/£>[/(*)] 
/£° [/(*)] 
f!f 2) [f(z)] 
A Pi) [f(z)] 



1 + 



k 2 + m' 



n 



rw + to - 1 n*) 



A*) 



1 / 2 W 



+ f /(*)/'(*) 









i) 








5) 



/' 2 (z) +1^-1) f\z) 



n 



11 



f%) + 



n 



1 / 2 (*) 



-/(*)/'(*) + '-f\* 

z z z 

fc 2 r 2 



(2.16) 
(2.17) 
(2.18) 



+ 'ih = 3 '--- D -i-(2.i9) 
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Integrating over k by using the formula 

-r s — r — - (2.20) 



k s d N k n+s-i^ N ~ 1)/2 ~ ( N + s-l\^{N + s 



VFT? T(N/2) 
we can see that 

Pi = -e, i = 3,...,D-l. (2.21) 

The v.e.v.'s ( |2.15| ) are divergent. To make them finite we introduce the cutoff function 7/^(7), 
which decreases with increasing 7 and satisfies the condition ip„ — ► 1, fj, — ► 0. To extract the 
divergent parts we will apply to the corresponding sums over I the summation formula [13] 

i f 00 K„(x) e - n7ri / (ix) - e n ™f{-ix) 
+7T / dx TTT A 2 • 2 - 22 

This formula is valid for functions f(z) satisfying the conditions 

\f{z)\ < e{x)e c \ y \ z = x + iy, c < 2, (2.23) 
f(z) = o(z 2 ^~ 1 ), z^O, (2.24) 

where e(x) — ► for x — ► 00. 

To evaluate the mode sum over Z in (fyD) in (gjp function /(z) we choose 

/(*) = z 3 f^[J n (zr/a)]^(z/a). (2.25) 

We will assume a class of cutoff functions for which ( |2.25| ) satisfies conditions ( [2.23D and (|2.24|) 
uniformly with respect to cutoff parameter ji (this is the case, for instance, for e -/n , fi > 0). 
Using the relation f$ [J n {-izr/a)] = e 2nni fi q) [J n (izr/a)] we see that the subintegrand of the 
first integral on the right hand side of formula (2.22) is proportional to ijj^iz/a) — ip^—iz/a). 
Consequently, after removing the cutoff — > 1) the contribution of the first integral will be 
zero. Hence, omitting this integral for the v.e.v. we obtain 

, = _i_ /> k + f (1 r 4 *HMW + (2 . 26) 

(27r) Ar + 1 a 3 i n f" oo \2i Vz 2 + k 2 a 2 + m 2 a 2 V 7 



K JaVk 2 +m 2 I n {z) V ' Z 2 - k 2 CL 2 - m 2 a 2 J ' 

where Xu(y) = bPuiw) + The second integral on the right of this formula vanishes 

in the limit a — > 00, whereas the first one does not depend on a. It follows from here that 

the latter corresponds to the spacetime without boundaries. This can be also seen directly by 

+00 

explicit summation over n using the formula J 2 ± p (z) = 1. For instance, in the case of the 

n=— 00 

energy density it follows from here that 

+oc ,s k 2 +m 2 

E ti e) [Uzr)] = l + !^F-. (2.27) 

z z 
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Using this relation for the term corresponding to the first integral on the right of ( |2.26j ) and 
removing the cutoff one has 



,(0) 



1 



2(2tt) jv + 1 
1 



^ dz z 3 fi £) [j n (zr)]^{z/a) 



2(2vr) 



N+l 



d N k 



n=—oo 
roo 



\lz 2 + k 2 



(2.28) 



dz z\J z 2 + k 2 + m 2 = — 



f d D - l \i /-= 



Hence summation formula ( |2.22| ) allows us to extract the contribution of the unbounded space 
without specifying the cutoff function. In the remaining part, the integration over k can be 
explicitly done using the formula 



k s g(z)dz 



7T 



N/2 



f ,jv, f°° k s g{z)dz it"/' n /N + sl. . , . ., 

(2.29) 

where B(x,y) is the Euler beta function (this formula can be proved by integrating over the 
angular part of d N )t and introducing polar coordinates in the plane (x = y/ z 2 — k 2 — m 2 , k)). 
Removing the cutoff for this part one obtains 

2 2 - D ir- D / 2 +^ 
a°T(D/2 - 1) 



(JV+a-l)/2 



QSUB = q-q 



(0) 



dz z 3 (z 2 



2 2 
a m 



D/2-2 KJz) 



In(z) 



F^[I n (zr/a)}, 



where I n (z) and K n (z) are the modified Bessel functions, and 



F^lf(z)} 

F^[f(z)] 
F^lf(z)] 



1 



m 



D 

1 
2 



1 



n 



+ l\f 2 ( z )-f\ z 



2£ 



/(*)/'(*) 



(2.30) 

(2.31) 
(2.32) 



2f-i 

s 2 



n 



D-l. 



IP n 2 
+ -±f(z)f'(z) - -f(z) (2.33) 

z Z A 

(2.34) 



In (|2.30|) g(°) are the components of the v.e.v. (OjvflXf |0m) over the vacuum |0m) corresponding 
to the unbounded D-dimensional Minkowskian spacetime. It can be seen that v.e.v. 's (|2~30D 
satisfy the continuity equation V^T? = 0, which, for the cylindrical geometry, takes the form 



dpi .1 , 
- r + -{Pi 
dr r 



P2) 



(2.35) 



For D = 3 formulae ( |2.30D coincide with the corresponding results in jjT(| for the spherical case. 

Formulae ( |2.30| ) can be obtained by another equivalent way, applying a certain second-order 
differential operator on the regularized Wightmann function and taking the coincidence limit. 
The calculation of the Wightmann function is of interest for its own sake, as this function 
determines the response of a particle detector moving through the vacuum under consideration 
@. For this reason here we derive formula for the regularized Wightmann function. To do this 
we apply summation formula (|2.22| ) to the sum over / in flOp . It can be shown that the first 
integral on the right of ( |2.22| ) will give the corresponding Wightmann function for D-dimensional 
Minkowskian spacetime, and the second integral vanishes. Integrating over k with the help of 
standard integrals, in the term corresponding to the third integral one obtains 



(Q\cp(xMx')\Q) 



(0 M \ip(xMx')\0 M ) - 



+00 in(4>-4>') f oo 



(2tt) d / 2 



dzz{z 2 -m 2 ) '*- 1 x (2.36) 
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fr ( \ J D 12-2 
x I n (zr)I n (zr') Kf -■ 



z* — 



l\2 



(r,, -r'A 2 -{t-t') 



1/2^ 



In(z) 



(r|| _ r j |)2 _ (t _ t0 2 



D/4-1 



It can be easily seen that for the case of a circle (D = 3) this formula coincides with the 
corresponding expression for the spherical case derived in fi~6|| . 

Let's consider vacuum densities (2.30) on the cylinder axis, r = 0. Now the only contributions 
come from the summands with n = 0, ±1 and one has 



l(r = 0) 



2 2- DlT -D/2 

a D T(D/2-l) 



3 / 2 2 2\ D / 2 ~ 2 f 

dz z lz — a m ) < u q o 



K (z) K x {z) 

+ U q \- 



Hz) 



Hz) 



(2.37) 



where 



u pi0 



1 - m 2 a 2 z 2 1 , 1 

h2f - -, it e i = 2£ - - 

L>-2 s 2' s 2 



(2.38) 



u P2 o 



^Pll — ^P2l 



The results of the corresponding numerical evaluation for the Dirichlet and Neumann minimally 
and conformally coupled scalars in D = 4 are presented in Fig|l|. 
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on the cylinder axis for minimally (left) and conformally (right) coupled Dirichlet (a, b) and 
Neumann (c, d) scalars in D = 4 as functions of ma. 



In the large mass limit, ma 3> 1, using the asymptotic formulae for the modified Bessel 
function for large values of argument in the leading order from ( |2.37| ) one obtains 



e (r = 0) 2px(r = 0) 



4£ - 1 



2 D-l^D/2-l a D 



(am) D / 2+1 e- 2am (25 B0 - 1) , ma > 1. (2.39) 



The v.e.v.'s ( |2.30| ) diverge on the cylinder surface (for a given n the z-integrals diverge as 
(a — r) l ~ D ', r — > a). The corresponding asymptotic behaviour can be found using the uniform 
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asymptotic expansions for the modified Bessel functions, and the leading terms have the form 



(D-m-QT(D/2) 

£ ~ ~P2 2 D- % D/2( Q _ r) D (2*B0-1) (2-40) 

(g-e c )r(D/2) 

Note that these terms are independent of the mass and Robin coefficients. Taking the limit 
a,r — * oo, a — r = const one obtains the leading terms for the asymptotic behaviour near the 



single plate considered in ref.[15|. The surface divergences in the renormalized expectation values 
for the vacuum EMT are well known in quantum field theory. In the case of D = 3 massless 
fields the corresponding asymptotic series near arbitrary smooth boundary are presented in 
refe.pl pa. 



3 Total Casimir energy inside a cylinder 

The unregularized integrated Casimir energy inside a cylindrical shell may be obtained by inte- 
grating ( |2.15 ) for q = e. Using standard formulae for the integrals involving the Bessel functions 
this energy can be presented in the form 



2 J \2ttJ Jxl l + k 2 a 2 + m 2 a 2 



where 

1 / L x N 



is the total Casimir energy inside a cylinder, evaluated as a sum of the zero-point energies for 
each normal mode of frequency. As we see, for the general Robin case, ([O]) differs from the 
total vacuum energy. The reason for this difference is the existence of an additional surface 
energy contribution to the volume energy, located on the boundary r = a — 0: 

To see this, note that there is a surface energy density contribution to the vacuum energy density 

T ( rf) = - (2^ - *(r - a + O^n^. (3.4) 
The corresponding v.e.v. can be evaluated using eigenfunctions ( |2.5| ): 

(OI^IO) = 5(r-a + 0)^0 £ 7t=== 2 Jn(V/^(V/°)- (3-5) 

The integration of this energy density leads to the expression ([O]) for the total interior surface 
energy. In (|3.5|) , applying to the sum over 2 formula ( 2.22 ), in analogy to Q2.30 ), the subtracted 
surface energy density can be presented in the form 

»» /sub 2 D -% D / 2 r(D/2-l)o D - 1 n £- M J mo V ; 

•-^W-i^k^V (3-6) 



i„ z V a/ \ a 
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This leads to the following subtracted surface energy 



E, 



(surf) \ 
in / 



(4f - 1) L u - 6 a 



D-3„2-D 



sub 2 d - 2 tt d / 2 - 1 T(D/2 



dz z' 



(3.7) 

We can find the subtracted total volume energy inside a cylinder integrating the energy density 



E 



(vol)\ 



L D-3 a 2-D 



TT E / dzz { 
J- ) „ J ma v 



2 2 2 

z — m a 



/sub 2 d ~ 2 tt d / 2 - 1 T(D/2-1) 
where, for a given function f(z), we have introduced the notation 



D/2-2 R n (z) 

40) 



(3.8) 



F v [/0)] = (4e-l>/0)/'0) 



z 2 — m 2 a 2 



f' 2 (z)-[l + ^\f 2 (z) 



n 



D-2 

The total interior vacuum energy can be obtained as a sum of the surface and volume parts: 

L D-3 a 2~D +oo «, 9 xD/2-l 



(3.9) 



(3 



ml SUB 



2 d - 1 tt d / 2 - 1 T(D/2) 



K n (z) 

In(z) 



n=—oo 
2 N 



E / M 



2 2 2 

z — m a 



(3.10) 



1 + 



n 



I 2 n(z) 



As we see, in the total vacuum energy the dependencies on the curvature coupling £ are can- 
celled out and we obtain ^-independent vacuum energy. We could have expected this result, 
as eigenfrequencies are independent of this parameter. The vacuum force acting from inside 
per unit surface of the cylinder can be found using the expression for the vacuum radial 

pressure: 



-^in — Pi \r=a- 



2 2- D7T -D/2 



a D T{D/2 



1 1 „— ™ J ma v 



2 2 A D/2-2K n (z) 



z —ma 



In(z) 



F^[I n {z)\, (3.11) 



with notation ( p. 32 ). For the massless scalar field and in the case of conformal coupling, £ = £ c 
the vacuum force and volume energy are related by 



E 



(vol)\ 



2vrL 



D-3 



I sub D — 2 



-a 2 R n . 



(3.12) 



This relation can be also obtained from the continuity equation ( 2.35| ) if we take into account 
the tracelessness condition for the corresponding energy-momentum tensor. 



4 Casimir densities outside a cylindrical shell 

Firstly, let's consider v.e.v.'s for the energy-momentum tensor in the region between two coaxial 
cylindrical layers with radii a and b, a < b. The corresponding boundary conditions are in form 

(A + B^\<p(x) = 0, r = a,b. (4.1) 

The corresponding eigenfunctions can be obtained from ( |2.5[) with the replacement 

Jniir) -> g n (ian r ) = J n (7r)Y n (ja) - Y n {~fr)J n (^a), (4.2) 
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where Y n (z) is the Neumann function. The functions chosen in this way satisfy the boundary 
condition on the inner cylinder. From the boundary condition on the shell r = b one obtains 
that the possible values of 7 are solutions to the equation 

Cf (77,7a) = J n ( 7 a)y n ( 7 6) - F„( 7 a)Jn(7&) = 0. (4.3) 

The corresponding roots will be denoted by 7a = cr nj /, I = 1, 2, 

The normalization coefficient is determined from condition ( [2.9D , where now the integration 
goes over the region between the cylindrical shells, and is equal to 

vr 2 7T n a6 (7«) 



Ha 



Auja{2-K 



iiV+1 ' 



Here we use the notation 



1 JKvz) 



-1 



A z + B\rj z z l - n 2 ) 



B 2 (z 2 



b 

a 



(4.4) 



(4.5) 



Using the corresponding eigenfunctions it can be seen that the v.e.v.'s for the energy- momentum 
tensor have the form ( 2,14| ), where the components are given by formulae 



7T- 



q{r) 



4(2 7 r) Ar + 1 a 3 



E 



a n i + k 2 a? + m 2 a 2 



/n 9) [5n(0-n,i,Crn,«r/a)], q = e,p U 



(4.6) 



with the functions fi q) [f{z)} defined as (|27T^)- (|27T9|) . 

To evaluate the sum over / we will apply the summation formula 



7T 2 °° f 



h(x)dx 



J„ 2 (x) + y n 2 (x) 

h{xe m l 2 ) + h(xe~ 



(4.7) 



K n (rjx) 



tlx 



lo K n {x) K n (x)I n (r]x) - K n {i]x)I n {x) ' 

where we have assumed that all zeros for the function C^ b (rj, z) are real. In case of existence of 
purely imaginary zeros we have to include additional residue terms on the left of formula ( f4.7|) . 
Now, by an evaluation similar to (2.26), one can see that 

dz 



(2tt) d - 

7T 
+ - 



z 3 fi q) [g n (z,zr/a)} 



+ 



(4.8) 



J 2 (z) + Y n 2 (z) V^ 2 + fc 2 a 2 + m 2 a 2 

z 3 fi q) [g n (iz,izr/a)} K n (rjx) / K n (x)dz \ 
4 JaVk 2 +m 2 Vz 2 - k 2 a 2 - m 2 a 2 K n (x)I n (r]x) - K n (rjx)I n (x) J ' 

To obtain v.e.v.'s for the exterior of a single cylindrical shell let's consider the limit b — > 00. It 
can be easily seen that in this limit the second integral on the right of formula (|4.8| ) tends to 
zero, whereas the first one is independent from b. It follows from here that the expressions 



1 



2(2vr) Z) - 2 a 3 



r poo 



dz 



z 3 fL 9) [g n (z,zr/a)] 



J 2 {z) + Y 2 {z) Vz 2 + k 2 a 2 + m 2 a 2 



(4.9) 



are v.e.v.'s for the EMT components for the exterior region of a single shell. To regularize this 
quantities we have to subtract the parts corresponding to the unbounded space. As we saw, the 
latter can be presented in the form ( 2.28 ). Using the identity 



fn q) [g n (z,zr/a)] _ ( } - -I V Jn ^ f^\H^(zr/a)] 

J 2 (z)+Y 2 (z) *« Vn{zr,a)\- 2 ^.^ i n [H n (zr/a)\ 



(4.10) 
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with Ht\z), a = 1,2 being the Hankel functions, one obtains 

„ = I f d N k y y r^^V/i u U) 

4(2vr) c - 2 a 3 7 n^^Jo Hn\z) \/z 2 + W + m¥" V ' 7 

Assuming that the function Hn\z), (Hn(z)) has no zeros for < argz < tt/2 (— vr/2 < 
argz < 0) we can rotate the integration contour for z by angle tt/2 for a = 1 and by angle —tt/2 
for (j = 2. The integrals over (0, ia^/k 2 + m 2 ) and (0, —iay/k 2 + m 2 ) cancel out. Introducing 
the Bessel modified functions and integrating over k with the help of formula ( p. 29 ) for the 
subtracted v.e.v. one obtains 

2 2-D^-D/2 +00 f00 ^ a a ^ D/2 _ 2 /n(z) 

QSUB 



where we use notations (2.31)-(2.34). As we see, these quantities can be obtained from the ones 
for interior region by the replacements I —* K, K I. As for the interior region the v.e.v. 
( [4.12| ) diverge at cylinder surface. The leading terms of these divergences are determined by the 
same formulae ( 2.4Cp with replacement a — r —* r — a. 

Let's consider the asymptotic behaviour of v.e.v (4.12) at large distances, r — > oo, for the 
massless case. Introducing a new integration variable y = zr/a, and expanding the subintegral 
in terms of r/a, we can see that the leading term of the asymptotic expansion comes from the 
summand with n = 0: 

qsuB = ^ ln !\ D1 / / dyy D ~ l F^[K Q {y% r ^ oo. (4.13) 



Y(D/2 - l)r D lnr/a J 

The integral in this expression can be evaluated using the formula for the integrals containing 
the square of the McDonald function (see, for instance, [{19(1). This leads to the following result 

D — 2 F 3 (D/2) £ — £ 

e~ (D - 2) P1 ~ -— „ ~ jD/2 ^ 1 t) -^-fa- , r -* oo. (4.14) 

For a conformally coupled scalar this leading term is zero and e ~ l/r-° +2 . 

Integrating the vacuum energy density over the region outside of a cylindrical shell we obtain 
the corresponding volume energy: 

P (vol)\ -L D -*a 2 - D ^ [°° ( 2 2 2A D/2-2 I n ( z ) 

Epvt > = n n — FTTn i — ; — ; r / / (iz2 z -mo - . , F v \K n (z)\, 

ext /SC/B 2 D ~ 2 7r^/ 2 - 1 r( J D/2 - 1) n^Jma V ' L WJ ' 

(4.15) 

where the functional F v [/(z)] is defined as fl3.9|) . The outside surface energy can be derived 
using the surface energy density given by the formula ( |3.4| ) with the replacement a — — > a + 
in the argument of the 5-function. This gives 

/ P (Buri)\ (4g - 1) f 30 , 2 / 2 2 2\ D / 2 " 2 4(f) „ 

(Ez-t ) = r n — n ,n i — ; — ; r > / dz z [z — m a ] -= — — K n (z) K„ (z) . 

(4-16) 

Now taking the sum of ( [4.15 ) and ( 4.16| ) for the total exterior vacuum energy one obtains 

L D-3 a 2-D +oo \X>/2-l 



jn(f) 



n=— oo 
,2" 



K'*(z)-\l + -f\K*{z) 



z 2 



10 



As we could have expected, the dependencies on the curvature coupling are cancelled. The 
expression for the radial projection of the vacuum force acting per unit surface of the cylinder 
from the outside directly follows from ( 4.12| ) with q = p\\ 



-2 2 - D ir- D / 2 



+00 



"ext 



~Pl \r=a+O z 



a D T(D/2 



IT 5: 



7 3/2 22 

dz z z — ma 



(4-18) 

where the function F? 1 is defined as ( 2.32| ). Now we turn to the case of a cylindrical shell with 
zero thickness, assuming that the coefficients in the Robin boundary conditions are the same 
for the interior and exterior regions. The total vacuum energies end force per unit surface can 
be obtained summing the corresponding quantities for these regions: 



( E )sUB - ( E '^)sUB + ( E ext)sUB J F - F in + F ext . 

Using the expressions for the interior and exterior quantities we have 



(4.19) 



v (AP 1) T D - 3 n 2 - D +oc 

ci(surf)\ _ Ifg ~ l ) L a 

/SUB 2 D - 2 n D l 2 - l T{D/2 - 1) J^^ 



dz z[z — ma' 



2 2\ D ' 2 - 2 



(l n (z)K n 
1-/3 V - . . - 



» 



zV n {z)K> n {z) 



(4.20) 



for the surface part of the vacuum energy, 

-L D - 3 a 2 ~ D 



(E) 



SUB 



2 d - 1 tt d / 2 - 1 T(D/2) n 



E 



dz 



< 2 2 2^/2-1 
[z — ma) 



2j3 + (z 2 + n 2 -p 2 



(i n ( z )k n {z) 
zi' n {z)k' n { z ) 



(4.21) 



for the total vacuum energy, and 
F 



_ 2 l-D n -D/2 

a°r(D/2 - 1) J< 



2 x D/2-2 



2/3 - 4£ + (z 2 + n 2 -(3 2 + 4£0) 



dz z (^z 2 — m 2 a 2 ^j 

(i n (z)k n { z ))' 



zI' n {z)K' n (z) 



(4.22) 



for the total vacuum force acting per unit surface of the shell. In these formulae we have 
introduced the notation 

f(z) = z' 3 f(z), P = A/B (4.23) 

for a given function f(z). The cases of Dirichlet and Neumann boundary conditions are obtained 
taking the limits (3 = 00 and (3 = 0. For the surface energies from ( 4.20 ) one obtains 

D-3„2-D +00 



E 



(surf)\( D ) 



I SUB 



p(surf)\( N ) _ (l-4g)I/W v 
/sub 2 d - 2 tt d / 2 - 1 T(D/2-1) ^ 

V ' / n= — cx 



, , 2 2 2\°/ 2 ~ 2 

dz z [z —ma 



For the analytical continuation of the expression on the right we note that 

+00 



E 



dr. z (^z 2 — m 2 a 2 ^j 



D/2-2 1 



B(D/2 - 1, 1 - D/2)(l + 2Ca(0)) = 0, 



(4.24) 
(4.25) 
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where denotes the Riemann zeta function (Cr{s) = /J n s ). As a result, we conclude that the 



n=l 



surface energy is zero for Dirichlet and Neumann scalars. Due to the relation ( 4.25| ) we can omit 
n- independent terms in the subintegrands of ( 4.20 ), ( [4.22[) without changing the values of the 
sums. We can use this to improve the convergence properties of the corresponding z-integrals. 

Let's present the total Casimir energy (4.21) in another equivalent form. For this, note that 
the functions f(z) = I n (z), K n (z) satisfy the equation 



z 2 f" (z) + (1 - 2f3)zf(z) - (z 2 + n 2 - (3 2 )f(z) 
Using this, it can be easily seen that 



0. 



2/3 + {z 2 + n 2 -(3 2 



(i n k n )' _ (r n K)' 



zV K' 



V K' 

n n 



+ 2(1-/3) 



(l n K n 



T K 



(4.26) 



(4.27) 



Thus, we obtain an equivalent representation for the total Casimir energy: 
—L 



(E) 



D-3 a 2-D 



SUB 



2 d - 1 tt d / 2 - 1 T(D/2) n 



E 



dz(z 



m 2 a 2 ) D l 2 ^ j-\n{l n {z)K n (z)) 



(4.28) 



This expression is divergent in the given form. A method to extract finite results will be explained 
below using the case of the D = 4 massless scalar field as an example. The corresponding results 
for arbitrary dimension and massive case will be presented elsewhere. 



5 Integrated Casimir energy per unit-length by mode sum 

Taking advantage of the usual prescription for evaluating the eigenfrequency sum (see. e.g. 
ref.|]]), the integrated Casimir energy per unit-length of an infinitely long cylinder is expressed 
as 

e c = ^Ccyl(CT) , (5.1) 

* <7->- 1 

where Ccyi(°") is the zeta function for the cylinder eigenfrequencies under the chosen boundary 
conditions. Note that, in order to simplify the procedure, we have not included the arbitrary 
mass scale typically used in these problems (in fact, it will not be necessary because the result 
will be finite). This zeta function can be related to the one for a circle, with the same boundary 
conditions and radius a, which we shall call Cdrc( s )> an d corresponds to the restriction of the ini- 
tial problem to a plane perpendicular to the cylinder axis. Because of the cylindrical symmetry, 
the eigenvalues of the wave operator — including mass — have the form (■f 2 + k 2 +m 2 ) 1 / 2 , where 
the 7's are the eigenfrequencies of the restricted problem and fceR. Then, in a spacetime of 
D = N + 2 + 1 dimensions, 



. , , f d N k / 2 ,2 2\- a / 2 1 r ( 2L 2^ 

Ccyl(-) = J J^W E (7 +k +m) = ^ p( _ §) 



Ccirc(<T - N) (5.2) 



with 



CcircOO = E(7 2 + ™ 2 r /2 = E d n E(7n, + m 2 )^ 2 , ln , = ^. (5.3) 

7 n I a 

In this section we are considering just one cylindrical surface, and these 7's denote now the 
zeros of J n (7a)-ffi 1 ^(7a) (the first factor coming from the interior propagation and the second 
from the exterior propagation). The n index corresponds to angular momentum, d n indicates 
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each degeneracy, and the I index labels the different frequencies for a given n. Expanding (5.3) 
around a = —1, one distinguishes two possibilities depending on N: 

1. Even N (N = 2p) 

Ccyi(a) = 22p+ 4 +1/2 r (-p - i) Cc irc (-2p - l). (5.4) 

2. Odd AT (JV = 2p+ 1) 

(-I)p 

Ccyl(CT) = 22p 7 rP+ 1 (p+l)! 
W< _2/ ' ~ 2 ' 1 (V(-l/2) - HP + 2))Ccirc(-2p - 2) + Ccirc'(-2p - 2) + 0{a + 1) 



CT+ 1 2 

"(5.5) 

Later on we will focus on D = 3 + 1 (i.e., p = case). Then, in situations where Ccirc( — 2) = 0, 
(|5.5|) reduces to 

Ccyl(s) = ^Ccirc'(-2). (5.6) 

Thus, one has to obtain the expansion of Ccirc(s) up to linear terms in {a + 1) = (s + 2) (in 
general, a + l = s + D — 2). By the Cauchy theorem, we may express the circle zeta function 
Ccirc(s) as the following contour integrations 



(5.7) 



Ccirc(s) = a s Y.dn^~. I dz(z 2 + a 2 m 2 )- s / 2 ^ In (j n (z)H n (1) (z) 
~ 2m Jc az V 

= a s J2 d n^-f dz{z 2 + a 2 m 2 )-^ 2 z\n(j n {z)H n ^{z) 
n 27TZ Jc v 

where C is a circuit in the complex z-plane enclosing all the zeros of J n {z)H r } l \z) (it is un- 
derstood that in this function we have now set a = 1). In this way, we are already including 
the contributions associated to internal and external propagation. Formulas (|5.7|) have the 
disadvantage of being valid for an s-domain in which Re s is positive, while we need to reach 
s = — (D — 2). Just in order to study the formal relations between these functions, we take 
the first form in ( |5.7|) and consider a specific circuit. Let C be the contour made of a vertical 
line along the imaginary axis avoiding — and leaving outside — the cut of the(z 2 + a 2 m 2 )~ s / 2 
function (placed between —iam and +iam), and a semicircle of infinite radius to its right. After 
realizing which parts yield no contribution, one finds 

Ccirc(s) = a s £ dn - sin (^) r dx (x 2 - a 2 m 2 )- s / 2 In (l n (x)K n (x)) . (5.8) 

IT \ Z J Jam <lX 



Now, using Q5.2p , one may formally write the integrated Casimir energy per lateral unit-length 

as 

- " - 2 »S-»r<r(%) ? d " £ dx <l2 " aVr "* I >n ( ' ;(I)A '" (I)) ' 

(5.9) 

Nevertheless, we have to stress that it is not valid to directly set s = — (D — 2) in order to obtain 
Eq/L d ~^. Analytic continuation to the left of the present s-domain is required first. In the 
example below, we explain a method for performing this type of process. Note, by the way, that 
if one 'naively' sets s = — (D — 2) in (|5.9|), expression ( 4.28j ) is recovered. 
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5.1 m — 0, D = 3 + 1 case 

In order to obtain the necessary analytic continuation to s = —2, we will follow a different path, 



which has been explained e.g. in ref. [21]. Its application will be illustrated for the case m = 0, 
D = 3 + 1 (therefore do = 1 and d n = 2 for n / 0). The mode sum is decomposed into sets of 
partial waves with a defined angular momentum n, which, for convenience, will be separately 
considered. Thus, we employ the notation 

oo 

CcircOO =E d nUs) = ]T ( n (s), (5.10) 
n n=—oo 

where, taking the second form in (5.7), 

oo „ 

Us) = = ^~ \ c dzZ ~ S ~ X 111 (jn(z)Hn (1 \z)) • (5.11) 

A first step in the extension of the s-domain to the left is the study of the asymptotic behaviour 
of the integrand for \z\ — > oo, given by J n (z)H n (z) ~ —iz/2 = / as (.z). In view of this, we 
subtract from and add to the integrand a term z~ s ln/ as (z). This move leaves us with 



f dzz-° In ( J n^)Hn\z) \ + j dzz „ s 



(5.12) 



Taking a contour C as described (now, since m = 0, the cut line has reduced to the point z = 0) 
we must conclude that the second integral vanishes, because its integrand has no singularity in 
the interior. As for the first integral, there are nonvanishing contributions only from the vertical 
parts of C. Parameterizing them adequately, we may write 



( n (s) = a s ^-sin(^\ J dxx s lnL n (x) 



2 

L n {x) = -- [I n (x)K n (x)\ 



(5.13) 



This formula is slightly better than the initial one, as its domain of validity has been shifted to 
the left: actually, it holds for — 1 < Res < 0. However, it is not good enough, as we have not 
arrived at s = — 2 yet. To this end, we shall have to perform even more subtractions. Before 
going on, we decompose Ccirc(s) into the contribution from nonzero angular momenta (n ^ 0) 
and the contribution from zero angular momentum (n = 0): 

t ci r C (s) = C£°(s) + C= (s). (5.14) 

When n / we may perform a rescaling x -> nx. Afterwards, summing the internal and 
external contributions, we find 

Cn(s) = a s n~ s -sin I — J J dxx~ s ~ x \\\L n (nx). (5.15) 

We shall actually assume B = 1 and take A as the only variable. This can be done whenever 
B 7^ 0. Therefore, the only case which cannot be described in this way is B = 0, A ^ 0, 
i.e., purely Dirichlet conditions, which have already been considered in ref.Q (formula (34) in 
that paper). Studying the uniform asymptotic expansions of the involved Bessel functions (for 
B 0) one realizes that 



LJnx) ~ 



(1 + x 2 ) 1 / 2 / \ 

\B 2 + next-to-leading termsj , for large nx. (5.16) 
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This behaviour motivates (for B = 1) the subtraction from and addition to the integrand in eq. 



( |5.15|) of a term x s 1 In 



;i + z 2 ) 1/2 



dxx s 1 In 



Using then the result 
f(l + x 2 ) 1/2 \ 



7T 



2ssin(^)' 



(5.17) 



one obtains 



Cn(s) = a s {n s -sin — 



7T 



7TS\ 



2 J 



dx x 



-s-l 



In 



(1 + x 2 ) 1 / 



l/2 L n(nx) 



+ — . (5.18) 



Hence, 



n=l 

a s < 2 — sin 
I * 

+Cfl(s)}. 



m=l 



dxx s 1 In 



(1 + x 2 ) 



IT^^n(ra) 



(5.19) 



All this comes from a formula valid for Res between —1 and 0. As already remarked, in order 
to work at s = —2, further subtractions need to be done to the integral in (|5.19| ). The result of 
applying such a technique is 



Ccirc ( s ) 



where 



sin 



7T 



TTS 
Y 



oo N/2 

2 ]T n~ s S n (N, s) + 2J2 Cr(s + 2j) ^W 2ij& B 

n=l j=X k 



S + k S 

2 '~2 



<5 n (iY, s) 



/•oo 




/ dxx~ s 


In 


Jo 





AT/2 



(1 + x 2 ) 1 / 



YH L n{nx) -2^2 



U 2j (t(x)) 



i=i 
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In formula (|5.20[) , the objects called U2j{t) are polynomials in t, where 

- 1 

t[X) ~ (1 + X 2 )V2' 

which are constructed step by step when obtaining the asymptotic approximation 



(5.20) 
(5.21) 

(5.22) 



In 



(l + s 2)l/2 



L n (nx , 



N/2 



U 2j {t{x)) 



n 



2.7 



(5.23) 



for large nx — in fact, this is like a further elaboration of ( |5,16D — . Then, the U2j,k' s are the 
present coefficients of these polynomials, in other words, 



^2j(t) = y^2 ? -,fc t k , 



(5.24) 



where the /c-index range corresponds to the existing coefficients. The whole process for obtaining 
this type of approximation has already been described in detail in refs.[^, ^TJ (see also jp2| ). 
By comparison, one may reason that these U2j(t)'s have to be the same as the even-j Uj' (t)'s 
in formula (3.37) of ref. [21], after the replacement a — > A. The N variable denotes the order of 
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the subtractions performed in order to remove negative-s poles from the initial integral. If N is 
large enough — and in our case this means just N > 2 — the S n (N, s) integral is finite at s = —2. 
When N is further increased, the value of the numerical integration S n (N,s) gets smaller and 
smaller while the algebraic n-sum (easier to evaluate) grows accordingly. In practice, we take 
N =8 or 10. 

The n = contribution has to be dealt with separately, as a step in the process for n / 
involved rescaling by n, and cannot be applied now. The same philosophy as in refs. 0, 1 will 
now be adopted. We will just perform the necessary subtraction to isolate the s = —2 divergence 
from the integral, leaving the rest to numerical evaluation. The result is expressed in the form 



3^00 = a' W ^ 

7T 



-R-o(s) - (a 2 



A + 3/8) i B (l+H, 



where the integral 



7*o 00 



J™ dxx~ s [ln|Z (a?)| + (A 2 - A + 3/8) t 2 (x) 



(5.25) 



(5.26) 



is, by construction, finite at s = —2. Therefore, TZq(— 2) can be calculated numerically for given 
values of A (B = 1). 

Thus, the circle zeta function is the sum of ( |5.20| ) and ( 5.25| ). At s = —2, ( 5.20|) has a 
nonzero finite value coming from the pole of B ^^p, — §) multiplied by the zero of ^ sin (^) at 
s = —2. This nonzero value depends also on the involved ^2j,2 coefficient (s) and on (r(s + 2j). 
Now, making a similar reasoning for ( 5.25| ), we realize that the s = —2 value of this expression 
is just — -^z{A 2 — A + 3/8). In fact, algebraic calculation shows that the finite value of ( 5.20| ) at 
s = —2 amounts to the same quantity with opposite sign. Therefore, Ccirc(— 2) = 0, and one can 
safely apply formula (|5.6| ). In these circumstances, the matter reduces to finding Ccirc' (—2) by 
expanding the sum of fl5,20| ) and (5.25) around s = —2. Once this has been done, the integrated 
Casimir energy per unit-length is simply 



(5.27) 



Results for a given ^4-range around A = (purely Neumann) are shown in Fig. || below, where 
we have plotted the numerical values of e c • a 2 as a function of A. 



6 Final remarks 

In the present paper, the vacuum — or zero-point — energy of scalar fields obeying Robin bound- 
ary conditions on cylindrical surfaces has been studied. Local quantities have been evaluated 
starting from vacuum expectation values of energy momentum tensors and applying generalized 
Abel-Plana summation methods supplemented with a suitable cutoff function. Integrated ener- 
gies seem to be more easily found from the eigenfrequency sum which determines the Casimir 
effect. In that case, zeta-function regularization has proven to be an adequate tool. However, 
both approaches constitute different aspects of one underlying idea. In fact, we have shown that 
the two methods lead to identical formal expressions — formulas ( 4.28j ) and (5J3) at s = —(D — 2). 



We have obtained an expression for the integrated energy inside the cylinder. In general, 
the total vacuum energy does not coincide with the integral of the volumic density, because 
there is a purely superficial contribution, located at the boundary itself, which has to be added 
to the volumic part. Moreover, the contributions associated to the modes propagating in the 
exterior of the cylinder has also been considered. In both cases, there is a part coming from the 
surface ('viewed' from inside, in the first, and from outside, in the second). The total Casimir 
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Figure 2: Total integrated Casimir energy per unit-length, multiplied by a 2 , for B = 1 as a 
function of the A coefficient, in a region around ^4 = 0. Note the presence of a minimum near 
A = —0.010. The value for A = (purely Neumann b.c.) is of —0.014176... in agreement with 
formula (44) of ref. f|] • 
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energy arises as the sum of all these contributions, and may also be obtained by eigenfrequency 
summation. As an example, we have chosen the massless case in D = 3 + 1 dimensions for 
application of the zeta function method. The result for a given parameter range is shown in fig. 
^, which displays the presence of a local minimum near (but not coinciding with) the point where 
Robin conditions become purely Neumann conditions. This indicates that, in a theoretical space 
of possible boundary conditions, some linear combinations are energetically privileged. 
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